Dynamics of the Measurement of Nuclear Spins in a Solid-State 

Quantum Computer 
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We study numerically the process of nuclear spin measurement in a solid-state quantum computer of 
the type proposed by Kane by modeling the quantum dynamics of two coupled nuclear spins on 31 P 
donors implanted in silicon. We estimate the minimum measurement time necessary for the reliable 
transfer of quantum information from the nuclear spin subsystem to the electronic subsystem. We 
also calculate the probability of error for typical values of the external noise. 



I. INTRODUCTION 

Recently, a new implementation of a,solid-state quan- 
tum computer was proposed by Kanea. Kane's idea is 
to realize the quantum qubits by using the nuclear spins 
of 31 P donors in silicon. These S=l/2 spins are known 
to exhibit very long relaxation times because of their 
weak coupling to the environment!]. At low tempera- 
tures, when there are no thermally excited electrons in 
the conduction band, this contact is mediated primarily 
by the magnetic field and the hyperfine interaction with 
s-electrons localized on the donor. However not all s- 
electrons are equally important for this mediation. The 
electrons in the inner shells are localized within one lat- 
tice cell and are also very weakly coupled to the environ- 
ment. Of the five outer (3p) electrons in a P atom (one 
more than Si), four form valence bonds with surrounding 
Si atoms, while the fifth is given off to the conduction 
band to form a loosely bound, hydrogen-like s-state in 
the field of the positively charged donor. Since the mo- 
tion of this electron is described in terms of the effective 
electron mass, which is small in semiconductors, its lo- 
calization radius can be as large as 30 A. Hence, its 
state may be controlled by relatively moderate electric 
fields. On the other hand, the exchange interaction be- 
tween loosely bound electrons located on different donors 
results in a weak indirect coupling between their nuclear 
spinstl, which coupling allows the implementation of, e.g., 
the quantum Control-Not operation. 

A schematic illustration of the quantum computer pro- 
posed by Kane is shown in Fig. [j] for the particular case 
of two-qubits. Two 31 P donors are implanted in silicon 
and subjected to an external dc magnetic field of about 
B = 2T. This creates a Zeeman splitting of the nuclear 
spin levels of about 3.5 x 10 7 Hz and a Zeeman splitting 
of electron levels of about 5.6 x 10 10 Hz (much smaller 
than the splitting between the ground and the lowest 
excited hydrogen-like states of electron, which isfm 15 
meV). The hyperfine coupling constant for Si: 31 P is 29 
MHz. Since this constant is proportional to the proba- 
bility of finding the electron near the nucleus, it can be 
decreased by effects that attract the electron away from 



the 31 P nucleus, which in the device depicted in Fig. |l| 
can be accomplished by applying a positive voltage to the 
"gates" labeled Ai and A2 . Similarly, the exchange inter- 
action between the electrons located on different donors 
may be tuned by applying positive or negative voltage 
to the J-gate in Fig. thereby changing the overlap of 
electron wave functions. 
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FIG. 1. Two P donors in silicon. The nuclear spins are 
coupled to the outer electrons by the hyperfine interactions, 
which can be controlled by the A gates. The electrons are 
mutually coupled via an exchange interaction, which can be 
controlled by the J gate. 

In this article, we do not consider the actual process of 
quantum computation but instead focus on the retrieval 
of the result after the computation has been completed. 
The weak coupling of the nuclear spins with their en- 
vironment, which is essential to avoid the decoherence 
that will spoil the quantum computation, makes this re- 
trieval a highly nontrivial task. In particular, it can- 
not be accomplished by means of current NMR methods, 
since their sensitivity is as yet insufficient for detecting 
the signal from a single nuclear spin. Recognizing this, 
Kane proposed a special measurement procedure based 
on transferring the information about the spin state of 
the system to its charge state, where the result can be 
measured. In the following sections, we simulate numeri- 
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cally the measurement procedure proposed by Kane and 
estimate its optimal duration and minimum probability 
of error. 



II. EIGENENERGIES AND EIGENSTATES 

The Hamiltonian of the full system is 

H = 2hbB{S Xz + S 2z ) - 2g n fi n B(h z + I 2z ) 

+4AiSiii+4^ 2 S 2 i 2 + 4JSiS 2 , (1) 

where Si and 1^ are electron and nuclear spin operators 
corresponding to donors i=(l,2), /is (/i„) is the Bohr 
(nuclear) magneton, g n — 1.13 is the nuclear g- factor, B 
is the external magnetic field, Ai is the hyperfine interac- 
tion constant for nucleus i, and J determines the strength 
of the exchange interaction between the electrons. 

The Hamiltonian H can be represented as a 16 x 16 
matrix in the basis of states with definite electron and 
nuclear spin projections. Since the total spin projection 
in the field direction is conserved, all possible states fall 
into five invariant subspaces, corresponding to its values 
-2, -1, 0, 1, and 2. In what follows, we will be interested 
only in the states with S z -£ + I z % = — 1, which are used 
for measuring of the nuclear spin states. Hence we may 
focus on a reduced basis of four states: | ||) e | |T)n, 
U)e| U)n, I IT)e| U)n, and | U) e | II)™- In this reduced 
basis, the Hamiltonian can be represented by the 4x4 
matrix 



H = 



J-2fi B B 2Ai 

J-2/i B B 

2Ai 2g nf i n B 

2A 2 23 




For simplicity, we henceforth assume that A\ — A 2 — 
A. Then the Hamiltonian is symmetric with respect to 
donors 1 and 2, and the eigenstates are either symmet- 
ric or antisymmetric with respect to interchanging them. 
The two symmetric states, \E±) and \E 2 ), have total elec- 
tron and nuclear spins Se = 1 and Is = 1. The corre- 
sponding eigenenergies are given by 



E x = g n ^ n B + J - n B B + yj (/ibB + g n fJ.nB) 2 + AA 2 , (3) 
E 2 = gnVnB + J — fi B B - y / (n B B + g n n n By + 4:A 2 , (4) 
and the corresponding (unnormalized) eigenvectors are 

\Ex) = -{2g n ^ n B + J -Si) jj) e | |t + Tl>n 

+2A\ |T + U)e| II>„, (5) 
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FIG. 2. Energy levels of the system of two coupled nuclear 
and electron spins versus exchange coupling J. The inset pro- 
vides an expanded scale view of the avoided crossing ("bot- 
tleneck") formed by E 3 (J) and i?4(J) curves at J = hbB/2. 

The Hamiltonian H also has two antisymmetric states, 
\E^) and I-E4), with total spin + Is = 1- The eigenen- 
ergies E3 and E4 are given by 



E z = g n iJ, n B -J- hbB + y/insB + g nf i n B - 2J) 2 + 4A 2 , (7) 
£4 = g„fi„B -J- v B B - yl {^bB + g nf i n B - 2J) 2 + 4A 2 , (8) 

and the corresponding (unnormalized) eigenvectors are 

\E 3 ) = -2A\ ||) e | |T - U>n 

-(2n B B- J + E 3 )\ IT - U>e|U>™, (9) 

\E 4 ) = -2A\ ||) e | |T - Tl>n 

~{2^ B B- j + e 4 )\ it - U)e| U>„, (10) 



Unlike the energies of symmetric states, £3 and E4 
are non-monotonic functions of J (see Fig. |J). If A 
were zero, the £3 (J) and £4 (J) curves would intersect 
at J = hbB/2. However, since these terms have the 
same symmetry, they repel to form an avoided crossing 
( "bottleneck" ) of width 4A at the would-be intersection 
point (see Fig. || inset). The corresponding eigenvec- 
tors also interchange character at this point. That is, 
as J increases from J < /is_B/2 to J > (XbB/2, \E$) 
transforms from | |T — TI)e| Ii)n with Se = and 
I s = 1 into I ||) e | |T - M) n with S s = 1 and I s = 0. 
In other words, the electron and neutron subsystems 
exchange their spins. Simultaneously, \E±) undergoes 
the inverse transformation from | H) e | IT — TI)n m to 
I IT — TI)e| II)n, with the opposite spin transfer. 



\E 2 ) 



H n B + J 
+2A\ |T + U)e 



-E 2 )\ ||) e 



IT + TI)r 



(6) 



The eigenenergies Ei and £ 2 are linear functions of J 
(see Fig. ^), and the structure of the corresponding eigen- 
states is independent of J (see Fig. ||). 
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FIG. 3. The components of \Ei) ,\Ei), |i?3), and \Ea) states 
versus J. 

One can thus distinguish between the singlet and 
triplet states of the nuclear subsystem as follows. Sup- 
pose that after the quantum computation has been per- 
formed with J = 0, the electron subsystem is in the | J.|) e 
state. Hence the whole system is either in \E 2 ) or in \Ei) 
state, depending on whether the nuclear subsystem is in 
the | |t + Tl)n or | |t - tl)n state. If the exchange pa- 
rameter is then adiabatically increased to J 3> pbB/2, 
the final electron subsystem is (respectively) a triplet or 
singlet state, thus allowing the information to be trans- 
ferred from the nuclear to the electron spin subsystem. 
To complete the measurement, one can distinguish be- 
tween these electron spin states through the difference in 
their charge properties, as described by KaneE. 



III. SIMULATIONS OF THE MEASUREMENT 
DYNAMICS 

To simulate the dynamics of measurement, we solve 
numerically the equation 



(ii) 



where p is the 4x4 density matrix of the system, 7i(t) 
is given by (^) with J(t) linearly increasing from to 
2/j.bB. We allow for a fluctuating correction 



5H = 45A 1 (t)S 1 I 1 + 46A 2 (t)S 2 I 2 



(12) 



to account for the voltage noise of the A-gates. The ran- 
dom classical quantities 5A% and SA 2 are assumed to have 
correlation time r c much shorter than all the dynamic 
time scales and zero averages. Following AbragamD, we 
rewrite the equation for the density matrix in the form 



dp _ i r 



^[Sn(t),[8H(t),p}], 



(13) 



where overlining denotes averaging over realizations of 
the fluctuations. We estimate the spectral density of 



fluctuations of A\ and A 2 using the spectral density 
of voltage fluctuations for room-temperature electronics 
S v = 1(T 18 V 2 /Hz and d(A/h)/dV ~ 30 MHz/V (see 
Rcf. |l|). The system is initially prepared in the \E^) state. 

Figure ^ shows the time dependences of occupation 
probabilities of the eigenstates \E\) - \Ei) during the 
measurement process. 
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FIG. 4. The time evolution of the occupancies of \Ei) - 
| £4} states during the measurement. The duration of mea- 
surement is T = 5.7 x 1CF 6 sec, the hyperfine interaction 

Hz, and the noise spectral den- 
Hz. 



constant is A/h = 2.9 x 10' 
sity is S A /h 2 = 3.5 x 1CT 3 



In this figure, we have taken the duration of measure- 
ment to be T = 5.7 x 1CP 6 sec, the hyperfine interac- 
tion constant to be A/h = 2.9 x 10 7 Hz, and the noise 
spectral density to be 5^//i 2 = 3.5 x 1CP 3 Hz. It is 
readily seen that increasing J at a finite rate results in 
a finite probability of exciting the system from \E^) into 
I-E3), creating one source of measurement error. Another 
source of error is noise-induced "escape" into the eigen- 
states I £4) and \E 2 ). Note that the system cannot be 
excited into these states in the absence of noise because 
of their different symmetry. Since the fluctuations of A-y 
and A 2 violate the symmetry with respect to donors 1 
and 2, they make this excitation possible. The different 
shapes of curves representing |Ci(i)| 2 and IC2 (*) | 2 are 
explained by the "restructuring" of the eigenstate |£Lt) 
state at J = pbB/2: the noise-induced transitions are 
accompanied by flipping one nuclear and one electron 
spin in opposite directions. The noise also contributes to 
the escape to \Es). 

Figure || depicts the dependence of final amplitude of 
^4) on the duration of measurement (solid curve). 
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FIG. 5. The amplitude of \E±) versus the duration of mea- 
surement for B = 2 Tesla and A/h = 2.9 x 10 7 Hz (solid 
curve). The dash-dotted and dashed curves show the same 
dependences for A/h = 2.9 x 10 8 Hz and A/h = 2.9 x 10 9 
Hz, respectively. 

For small T, | C4 1 2 — > because the short-time dy- 
namics of the electrons is only weakly affected by the 
hyperfine interaction, and the electron subsystem retains 
its spin during the course of the measurement. For suf- 
ficiently large T, P4 tends to 1 according to the law 
1 — \C±\ 2 oc exp(— T/t). For representative values of 
B = 2 Tesla and A/h = 2.9 x 10 7 Hz, one obtains 
t- 1 = 5.7 x 10 5 sec" 1 . 

To understand why r _1 is three orders of magnitude 
smaller than A, we calculated the same curve for A/h = 
2.9 x 10 8 Hz and A/h = 2.9 x 10 9 Hz (dash-dotted and 
dashed curves). Figure |(] shows the logarithmic plot of 
t _1 vs A. 



from noting that T/t ~ AE/hAt, where AE oc A and 
At, which is the time of passage of the system through 
the "bottleneck" region, is also proportional to A, since 
J is ramped up linearly in our simulation. 

If it were not the gate-voltage noise, the measurement 
error could be made vanishingly small just by increasing 
the duration of measurement. However, because of noise, 
the error passes through a minimum as a function of du- 
ration and then increases again. The details of behavior 
of the measurement error are shown in Fig. (7| 




FIG. 7. The measurement error, 1 — | C4 1 2 , versus the dura- 
tion of measurement for different levels of gate- voltage noise: 
S A /h 2 = 3.5 x 1(T 3 Hz (*), 3.5 x 1CT 4 Hz (solid line), and 
3.5 x 1(T 5 Hz (o). 



The minimum error (1 — \C±\ 



is proportional to 



the noise spectral density Sa (see Fig. ||) and is of the 
order of 1CP 6 for typical values of noise. 
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FIG. 6. The characteristic time r versus the hyperfine in- 
teraction constant A. 



By scaling A, we found that r _1 oc A 2 /2^BBh. This 
result can be readily explained. The width of the gap 
between the states \E±) and I.E3) is AE = A A (see Fig. 
2, inset), and the corresponding interval of J is also of 
the order of A. Hence the result that t _1 oc A 2 follows 
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FIG. 8. The minimum attainable measurement error ver- 
sus the noise spectral density. 
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